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as ■ ABSTRACT 

^. 

The Hawaii Imaging Fabry-Perot Interferometer (HIFI) on the University of Hawaii 
""T I 2.2m telescope was used to map the Ha and [O III] A5007 emission-line profiles across 

pLn ' the entire disk of the edge-on Sb galaxy NGC 4388. We confirm a rich complex of 

l/^ ' highly ionized gas that extends ^ 4 kpc above the disk of this galaxy. Low-ionization 

gas associated with star formation is also present in the disk. Evidence for bar 
streaming is detected in the disk component and is discussed in a companion paper 
(Veilleux, Bland-Hawthorn, & Cecil 1999; hereafter VBC). Non-rotational blueshifted 
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\^ , velocities of 50 — 250 km s are measured in the extraplanar gas north-east of the 

m 

(^ ' A redshifted cloud is also detected 2 kpc south-west of the nucleus. The velocity 

^^ ■ field of the extraplanar gas of NGC 4388 appears to be unaffected by the inferred 

■^ ■ supersonic (Mach number M ~ 3) motion of this galaxy through the ICM of the Virgo 

o . cluster. We argue that this is because the galaxy and the high-|z| gas lie behind a 

Q , Mach cone with opening angle -^ 80° . The shocked ICM that flows near the galaxy has 

|H I a velocity of ~ 500 km s~^ and exerts insufficient ram pressure on the extraplanar gas 

~^ \ to perturb its kinematics. We consider several explanations of the velocity field of the 

extraplanar gas. Velocities, especially blueshifted velocities on the N side of the galaxy, 
vj are best explained as a bipolar outflow which is tilted by >12° from the normal to the 

j_j ■ disk. The observed offset between the extraplanar gas and the radio structure may be 

due to buoyancy or refractive bending by density gradients in the halo gas. Velocity 
substructure in the outflowing gas also suggests an interaction with ambient halo gas. 
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1. Introduction 

The edge-on Sb galaxy NGC 4388 is projected <1° from the core of the Virgo cluster 
(Biggnelli, Tammann, & Sandage 1987; Yasuda, Fukugita, & Okamura 1997), but its systemic 
velocity of ~ 2,525 km s^^ is ~ 1,500 km s^^ above that of the Virgo cluster (e.g., Rubin, Kenney, 
& Young 1997; Yasuda et al. 1997). Nevertheless, studies by Chamaraux et al. (1980), Giraud 
(1986), Kenney & Young (1986, 1989), and Cayatte et al. (1990, 1994) show that the disk of 
NGC 4388 is deficient in HI; it appears truncated at ~ 10 kpc, well inside the optical radius. 
This suggests that NGC 4388 is experiencing the effects of ram-pressure stripping by the densest 
portions of the intracluster medium (ICM) in the Virgo cluster. Distance estimates using the 
Tully-Fisher relation suggest that NGC 4388 lies close to the mean distance of the M87 subcluster 
in the Virgo cluster, i.e. near the core of the Virgo cluster (Mould, Aaronson, & Huchra 1980; 
Sandage & Tammann 1984; Pierce & Tully 1988; Yasuda et al. 1997). In the rest of this paper, 
we will assume that NGC 4388 is indeed a member of the Virgo cluster and will adopt a distance 
of 16.7 Mpc (Yasuda et al. 1997), yielding a scale of 81 pc/". 

NGC 4388 was the first Seyfert galaxy discovered in the Virgo cluster (Phillips & Malin 1982). 
Nuclear activity has been detected at nearly all wavelengths. The morphology of the nuclear radio 
source suggests a collimated AGN-driven outflow (Stone et al. 1988; Hummel &: Saikia 1991; 
Kukula et al. 1995; Falcke, Wilson, & Simpson 1998). The nucleus hosts a strong hard X-ray 
source with power-law photon index T = 1.5, that is heavily absorbed by a column density A'^^ 
= 4 X 10^'^ cm^^ (Hanson et al. 1990; Takano &: Koyama 1991; Lebrun et al. 1992; Iwasawa et 
al. 1997). The large equivalent width of the iron Ka line in the ASCA spectrum of the nucleus 
(Iwasawa et al. 1997) and the detection of broad recombination-line emission outside (Shields & 
Filippenko 1988, 1996) but not at the nucleus (e.g., Miller 1992; Goodrich, Veilleux, & Hill 1994; 
Veilleux, Goodrich, & Hill 1997) are both consistent with a powerful hidden AGN in NGC 4388. 

At optical wavelengths, NGC 4388 has been known for some time to present extended line 
emission (e.g.. Ford et al. 1971; Sandage 1978; Cohna et al. 1987; Pogge 1988; Corbin, Baldwin, & 
Wilson 1988). A rich complex of ionized gas extends both along the disk of the galaxy and up to 
50" (4 kpc) above that plane. The extraplanar gas component has considerably higher ionization 
than the disk gas, and appears roughly distributed into two opposed radiation cones that emanate 
from the nucleus (Pogge 1988; Falcke et al. 1998). Detailed spectroscopic studies strongly suggest 
that this component is mainly ionized by photons from the nuclear continuum (Pogge 1988; Colina 
1992; Petitjean & Duret 1993). 

Several origins of this extraplanar gas component have been suggested: (1) it is tidal debris 



from a past interaction with another galaxy (Pogge 1988); (2) supernovae have ejected it into the 
halo of the galaxy (e.g., Corbin et al. 1988); (3) ram pressure from the ICM has elevated some 
of it from the disk (Petit- Jean & Durret 1993); (4) it has been ejected from the nucleus by the 
mechanism that is powering the radio jet (e.g., Corbin et al. 1988). 

This ambiguity stems from the weak constraints on the kinematics of the extraplanar gas 
that arise from the poor spatial coverage of the spectra used in previous studies. In this paper we 
use a Fabry-Perot spectrometer to provide complete two-dimensional coverage of the kinematics 
of the line-emitting region in NGC 4388. After a brief description of the data analysis methods 
(§2), we present the distribution and velocity field of the Ha and [O III] A5007 line-emitting gas 
derived from the data (§3). These data allow us to split the line-emitting gas into two distinct 
kinematic components, one associated with the galaxy disk and the other associated with the 
extraplanar gas. In §4 we derive dynamical parameters of the extraplanar gas and use them to 
constrain the origin of this component. The important issue of ram-pressure stripping by the ICM 
is addressed in the context of our results. In this section, we sketch a picture of NGC 4388 that 
accounts for the effects of both the supersonic ICM and nuclear activity (Fig. 5). Our conclusions 
are summarized in S5. 



2. Observations and Data Reduction 

The Hawaii Imaging Fabry-Perot Interferometer (HIFI; Bland & Tully 1989) was attached at 
the Cassegrain focus of the University of Hawaii 2.2 m telescope at the Mauna Kea Observatory to 
obtain in photometric conditions 54 seven-minute exposures of redshifted Ha and 38 eleven-minute 
exposures of redshifted [O III] A5007. The Ha spectra were taken on 1990 May 31 and 1990 
June 1 with a thick epitaxial GEC CCD (5 e^ rms readout), resulting in a pixel scale of 0'.'66 
pixel^^. A finesse 70 etalon with free spectral range 57 A produced a velocity resolution of 40 
km s~^. The [O III] spectra were acquired on 1992 March 12 with a blue-sensitive 1024 x 1024 
Tektronix CCD and a finesse 60 etalon having a free spectral range of 61 A. The spatial scale and 
spectral resolution of these data are 0'.'85 pixel"^ and 60 km s~^, respectively. Order separating 
filters with a flat-topped transmission profile centered at 6614 A (55 A FWHM) and 5030 A (36 
A FWHM) passed only one etalon order so true emission line profiles of Ha and [O III] A5007 
could be synthesized. Wavelength calibration exposures were made of diffuse Neon (Ha cube) and 
Mercury-Cadmium ([O HI] cube) sources. A current-stabilized, flat continuum source served to 
flatfield the images. Finally, the photometric standard stars HD 94808, rj Hya, and 6 Vir were 
observed through the same setup as the object to provide an absolute flux calibration. 

These data were reduced using the methods described in Bland & Tully (1989), Veilleux et al. 
(1994), and Bland-Hawthorn (1995). However, the Ha data suffered from several cosmetic defects 
associated with the aging GEC CCD. While only a fixed bias offset needed to be subtracted from 
the [O HI] data, row-by-row bias subtraction was done on each Ha frame to remove vertical 
striations of the bias pattern. After interpolating over the columns with poor charge transfer, the 



images of each night were stacked in a cube in increasing order of wavelength. Each data set was 
flatfielded using a similarly constructed whitelight cube. This method produced good flat fields 
for the [O III] spectra. Large-scale illumination gradients in the sky background that persisted in 
the corrected Ha cube were removed by dividing the spectra by a spatially smoothed sky frame. 
This frame was obtained by summing all of the frames in the data while interpolating through 
regions with Ha emission from the galaxy. The photometric accuracy resulting from this method 
is discussed in §3.1. The rest of the reductions was the same for both Ha and [O HI] cubes: spatial 
registration of the frames, cosmic ray removal, correction for air mass variations, sky continuum 
and line subtraction, correction for the curvature of the isovelocity surfaces ("phase calibration"), 
and flux calibration. The OH sky line at 6604.13 A (Osterbrock & Martel 1992) was also used for 
the velocity calibration of the Ha data. 

For the sake of simplicity, the observed emission-line profiles were approximated by simple 
Gaussian functions (regions where the emission-line profiles differ from a simple Gaussian will be 
discussed in §3.2). The best fits were determined by the least-chi-squared method on spectrally 
smoothed emission line profiles using a 1/4 - 1/2 - 1/4 spectral filter (Hanning smoothing). 
Gaussian smoothing of each image with a = 1.0 pixel improved the sensitivity to fainter features. 
In the end, the centroid, intensity, and width were successfully derived from ~ 5,000 Ha spectra 
and ~ 7,000 [O HI] A5007 spectra. 



3. Empirical Results 

3.1. Distribution and Excitation Properties of the Ionized Gas 

Figure 1 presents the distributions of Ha and [O HI] A5007 line emission that we derived 
from the Gaussian fits; the underlying starlight has been removed. The 1-a detection limit of 
our data corresponds to a surface brightness of ~ 7 and 3 x 10~^^ erg cm~^ s~^ arcsec"^ for 
Ha and [O HI] A5007, respectively. The fiux calibration of our data cubes was checked against 
published long-slit spectrophotometry (e.g., Pogge 1988; and Petitjean & Durret 1993). The Ha 
and [O HI] fiuxes agree with published ones to within --^ 25 - 50%. Note that the Ha and [O III] 
fluxes were not corrected for Galactic reddening because extinction is negligible toward NGC 4388 
[E(B-V) = 0.028; Burstein &; Heiles 1978, 1984]. The patchy dust lanes running roughly east- west 
in the optical continuum images of this galaxy (e.g., Phillips &: Malin 1982) suggest highly 
variable internal extinction which may affect the apparent distribution of the ionized gas near 
the galactic plane. No attempt was made to correct for this effect. Dilution of the Ha fluxes 
by underlying Balmer absorption is negligible in NGC 4388. Moreover, the Balmer absorption 
profile is dominated by pressure broadening in the atmospheres of early-type stars, making this 
feature indistinguishable from a slowly varying continuum over the narrow spectral bandpass of 
our instrument. 

The most prominent features in the Ha image are the strings of H II regions in the shape of 



two 'spiral arms' within the disk of the galaxy. Strong Ha disk emission is detected out to --^ 5 kpc 
east of the nucleus (its extent is limited by the velocity coverage of our data cube). The bright 
H II region complexes in this part of the disk each have Lhq- ~ IQ^^'^ - 10^^'^ erg s~^. Emission 
near the detection limit of our Ha data is observed projected north of the plane of the disk along 
PA = 20° - 50° and 1-4 kpc from the nucleus. In contrast, the [O HI] map is dominated by 
nuclear line emission. Emission from the spiral arms and from gas above the galaxy disk is also 
visible but considerably fainter. Its large projected |z| distance implies that the gas is well above 
the plane of this highly inclined (|i| = 78° it 1°; cf. §3.2) galaxy. If this gas was in the disk, it 
would lie ~ 20 kpc from the nucleus which greatly exceeds the visible radius of the galaxy {R25 = 
14 kpc from de Vaucouleurs et al. 1991). 

The differences between the Ha and [O HI] line emission distributions is particularly evident 
in Figure 2, which maps the [O HI] /Ha flux ratios. This map fundamentally compares the 
ionization level of the line-emitting gas. Energies of 35.5 eV and 54.9 eV are needed to produce 
0^+ and 0+++, i.e. significantly larger than the 13.6 eV needed to ionize hydrogen. Uncertainties 
on ratios range from 25% in the bright H II regions in the disk to >50% in the fainter high-|z| 
gas (recall that the extraplanar Ha emission is near the detection limit of the data). The 
low-ionization gas is found in the bright H II regions of the spiral arms. These H II regions have 
[O HI] A5007/Ha ratios of ~0.1, typical of star- forming regions with low ionization temperature 
(~ 40,000 K; Evans & Dopita 1985) or high metal abundance (0/H ~ 3 x O/H©; McCah, Rybski, 
& Shields 1985). The high-ionization gas concentrates into two regions north and south of the 
nucleus. The southern extension has a conical shape with apex near the nucleus. Typical values 
for the [O HI] A5007/Ha ratio in the nucleus and in the ionization cones are 1 to 5, consistent with 
photoionization by the AGN (Pogge 1988; Colina 1992; Petitjean & Durret 1993). This ionization 
map is qualitatively similar to previously published [O HI] A5007/(Ha + [N II] AA6548, 6583) 
ratio maps obtained with narrow-band interference filters (Pogge 1988; Corbin et al. 1988; Wilson 
& Tsvetanov 1994; Mulchaey, Wilson, &: Tsvetanov 1996; Falcke et al. 1998). Note, however, that 
the present line ratio map was derived directly from our spectra and is therefore not subject to 
the large uncertainties in the line separation and continuum subtraction which often affect line 
ratio maps derived from narrow-band images. 



3.2. Kinematics of the Extraplanar Material 

Figure 3 presents the barycentric velocity fields that are derived from the Ha and [O HI] 
A5007 spectra; uncertainties range from ~20 km s~^ in the bright line-emitting regions to ^100 
km s~^ in the fainter areas. The overall velocity field of the [O III]-emitting gas resembles that 
of the Ha-emitting gas. In both, the velocity field in the disk is characterized by a large-scale 
east-west gradient indicative of rotation in the plane of the disk. The rotation curve derived from 
our data is highly consistent with earlier results (cf. Rubin, Kenney, &: Young 1997 and references 
therein). Non-Gaussian Ha and [O HI] A5007 profiles are found near the nucleus (R ^ 10 "), again 



confirming earlier results (lye & Ulrich 1986; Colina et al. 1987; Ayani & lye 1989; Veilleux 1991; 
Rubin et al. 1997). 

Figure 3 shows a kinematic dichotomy between the disk gas and the extraplanar gas. To 
more objectively identify which portion of the line-emitting gas does not partake in disk rotation, 
we have constructed kinematic models to simulate the observed velocity field in the disk. These 
models are described in more detail in VBC. Because the focus of the present paper is on the origin 
of the extraplanar gas in NGC 4388, only the main results of our analysis of the disk velocity field 
will be discussed here. 

To minimize possible effects of the active nucleus, our analysis of the disk component focused 
on the Ha velocity field. To summarize our results. Figure 4 compares the observed Ha velocity 
field with our best-fitting axisymmetric disk model and our best-fitting disk model that includes 
bar streaming motions. Compared to the model with purely circular motions, the velocity residuals 
from the model incorporting elliptical streaming are significantly smaller throughout the disk. The 
residuals have a symmetric Gaussian distribution with dispersion of only 10 km s~^. This is an 
excellent fit considering that hydrodynamic effects (e.g., strong shocks along the bar; Athanassoula 
1992; Finer, Stone, & Teuben 1995) are ignored in our simple kinematic model. The parameters 
of the model with elliptical streaming are listed in Table 1 and discussed in more detail in VBC. 
The strong bar streaming motions prevented us from deriving a rotation curve directly from the 
data; instead, we used a smooth parametric function to model the rotation curve (cf. eqn. 3 in 
VBC). Note that the inclination of the disk, i = —78°, implies that the north rim of the disk is 
the near side (i.e. we are looking at the galaxy from under the disk). 

Our kinematic model suggests that most of the gas in the disk of NGC 4388 partakes in 
bar-forced orbital motion around the nucleus. There are, however, two notable exceptions. The 
Ha and [O III]-emitting gas with R ^ 1.2 kpc and FA = 10° - 70° has velocities which appear 
blueshifted by 75 - 250 km s^^ relative to systemic velocity. The location and velocity field of this 
gas implies that it is elevated above the disk of the galaxy. 

No obvious kinematic trends with height above the disk or with galactocentric radius are 
seen in the NE complex. However, a correlation appears to exist between the velocity field and 
the surface brightness of the gas in the NE complex. The most extreme velocities in this gas 
complex appear in the bright cloud located near the nucleus at R ~ 2 kpc and FA ~ 55° (v = 
- 250 km s~^; the 'NE cloud' in the nomenclature of Fogge 1988). Double-peaked [O HI] profiles 
with velocity splitting ~140 km s^^ are seen near the north-east edge of that cloud. The much 
fainter gas directly east of this NE cloud has velocities which average 175 km s""*^ larger (-75 
km s~^ relative to systemic) than that of the cloud. The trend between surface brightness and 
blueshifted velocities is also seen in the cloud complex at R ~ 3.3 kpc and FA ~ 35°. Here, the 
brighter clouds have velocities which are ~ 60 - 100 km s~^ more blueshifted than the diffuse gas 
surrounding them. 

The bright [O HI] cloud south-west of the nucleus at R ~ 1.8 kpc and FA = 215° also 



has anomalous velocities relative to the disk. This cloud is not detected at Ha, implying 
[O III] A5007/Ha ^ 5. Substructures (shoulders, line splitting) of order ~ 200 km s^^ are seen in 
the emission-line profiles on the north-east side of that cloud (also detected by Corbin et al. 1988). 
A comparison with the modeled disk velocity field indicates that this cloud has an average redshift 
of ~ 100 - 150 km s^^ relative to the underlying disk gas. Such supersonic motion through the 
disk ISM is highly unlikely; therefore, this cloud is probably in front of the disk. Its radial velocity 
relative to systemic is only --^20-30 km s^^. 

Contrary to Petitjean & Durret (1993), we see no general trend for gas at large radii to have 
negative velocities relative to the galaxy rotation. The probable origin of this discrepancy is the use 
of different kinematic models to reproduce the observed velocity field. As we demonstrated in §3.1, 
elliptical streaming is important in the disk of NGC 4388 but this effect was not incorporated in 
the axisymmetric models of Petitjean & Durret (1993). We feel that the complete two-dimensional 
coverage of our spectra and our use of slightly more sophisticated models put our conclusions on a 
more solid basis. 



4. Discussion 

4.1. Derived Properties of the Extraplanar Material 

A clue to the origin of the extraplanar gas is provided by the energetics of this material. In 
this section, we estimate the dynamical time scale, mass, and kinetic energy of the extraplanar gas 
in NGC 4388. 

An upper limit on the dynamical time scale of the extraplanar gas can be estimated from its 
observed velocity and hnear extent: tdyn ~ R/V ^ 2 x lO'^ (-R/4 kpc) (yobs/200 km s^"*^)^^ yr. 
The ionized mass above the disk follows from the [O III] A5007 flux. Other researchers have 
estimated the density and temperature of the extraplanar gas. Pogge used the [S II] A6731/A6716 
flux ratio to find densities in the range 100 - 1,000 cm~^. However, densities <50 cm~^ are 
claimed by Colina (1992) and Petitjean & Durret (1993) for the NE complex. Estimates for the 
temperature range from ~ 10,000 K (Petitjean & Durret 1993) to 15,000 K (Cohna 1992). In 
the following discussion, we parametrize the mass in terms of the density and use an electron 
temperature of 10^ K. 

The total [O III] intensity of •^ 2.5 x 10^^'^ erg s^^ cm^^ beyond R = 1.2 kpc and within PA 
= 10 - 70° yields an ionized mass in the NE complex of M(NE) Ri 4 x lO'' X^^ N;^2 ^0' where X 
is the fraction of doubly ionized oxygen and Ne,2 is the electron density in units of 100 cm^'^. In 
this calculation, we used a solar mass fraction for oxygen and a five-level atom approximation to 
estimate the emission coefficient (McCall 1984). A lower limit on the bulk kinetic energy of the 
north-east gas, £'buik(NE), may be derived by summing the bulk kinetic energy over each pixel in 
this region (Sj 1/2 5mi vf, where 6mi and Vi are the mass and observed radial velocity at pixel 



i). We find £'buik(NE) ^ 6 x 10^^ X^^ N^2 ^^gs. The "turbulent" kinematic energy of this gas, 
-E'turb(NE), may be derived from Sj 5mi af, where ai = Vfwhm/1-67; we get £'turb(NE) ?« 1 x 
10^^ X^^ Ng"2 ergs. Following a similar procedure for the SW cloud, we find M(SW) ?a 1,000 X^^ 
N~^ Mq, Ebuik(SW) ;^ 6 X 10^9 X-i N"^ ergs and Eturb(SW) « 3 x 10^° X'^ N"^ ergs. Using X 
= 1 and Ne,2 = 0.1, the total mass and kinetic energy of the extraplanar gas are thus ~ 4 x 10^ 
Mq and -Ekin = -E'buik + -E'turb ^ 1 ^ 10^^ crgs. Using our estimate of the dynamical timescale, we 
derive a lower limit on the time-averaged rate of kinetic energy injected into the extraplanar gas: 
dE]iin/dt ~ -Ekin/idyn ^ 2 x 10"^^ N~2 erg s~^. We must emphasize that the energetics derived 
in this section only take into account the optical component of the extraplanar gas. Entrainment 
of the molecular gas from the disk could significantly increase the energetics of the extraplanar 
material. 



4.2. Origin of the Extraplanar Material 

Our results can constrain the origin of the extraplanar gas in NGC 4388. In this section, the 
predictions of each of the four scenarios introduced in §1 are discussed and compared with our 
data. 



4-2.1. Tidal Debris from a Past Galaxy Interaction 

In analogy with the Magellanic Stream (Mathewson, Schwarz, and Murray 1977; Lin, Jones, 
& Klemola 1995), the extraplanar material in NGC 4388 may represent the tidal debris from a 
recent encounter with a small gas-rich galaxy (Pogge 1988). In this scenario, gas tidally stripped 
from the dwarf galaxy appears when ionized by the radiation cone of the active nucleus. This 
scenario has been suggested by Tsvetanov &: Walsh (1992) to explain the narrow ionized "arcs" 
in NGC 5252. In this SO galaxy with an AGN, the ionized arcs are invisible outside sharp-edged 
conical sectors. However, the arcs have kinematics (Morse et al. 1998) and space distribution that 
continue into H I gas (Prieto & Freudling 1996) on a ring inclined to the SO disk that is biconically 
illuminated, hence ionized, by nuclear photons. 

Numerical stellar dynamical simulations of the accretion of a dwarf galaxy by a large disk 
galaxy show that the satellite sinks in only ~ 1 Gyr and causes considerable damage to the disk's 
inner regions (Walker, Mihos, & Hernquist 1996). Satellites on bound polar orbits experience 
severe orbital decay, settling into an orbit coplanar with the disk after only one or two orbits 
(Walker et al. 1996). The well-behaved velocity field of the disk in NGC 4388 and the presence 
of gas well above the disk therefore tightly constrains the epoch of the purported encounter, 
i.e. ^ 10^ yrs. In this scenario, the velocity field of the tidal debris is expected to reflect the 
orbital motion of the original dwarf galaxy, showing trends with radial distance from the nucleus, 
the center of the potential dwell. These trends are absent in our data. Moreover, it is difficult 



to explain in this scenario the velocity - surface brightness correlation and the steep velocity 
gradients observed in the NE complex (§3.2). 



4-2.2. Interaction with the Intracluster Medium 

NGC 4388 is one of the most Hl-deficient spiral galaxies in the Virgo cluster (Chamaraux 
et al. 1980; Cayatte et al. 1990; 1994). A group III galaxy in the nomenclature of Cayatte et 
al. (1994), it has a severely truncated HI disk with an Hl-to-optical size ratio -Dni/^opt = 0.7. 
NGC 4388 is located very close to the center of the Virgo cluster. Taking the position of the core 
to be the peak position of the luminosity density, a = 12h 25.3m and 6 = +13° 18' and adopting 
16.5 Mpc for the core distance (Binggeli, Tammann, &: Sandage 1987; Yasuda et al. 1997), we find 
that NGC 4388 is <300 kpc from the center of the Virgo cluster. Using Vvirgo ~ 1,000 km s~^ for 
the recession velocity of the cluster (Yasuda et 1997), NGC 4388 is plunging edgewise at -1,500 
km s~^ into the core. It is therefore a prime candidate for ram pressure and turbulent viscous 
stripping by the ICM. 

Ram-pressure stripping has been proposed by Petitjean & Durret (1993) to explain the 
velocity field of the extraplanar material in NGC 4388. In this scenario, the line-emitting gas 
in the NE complex may represent (1) gas that has cooled and condensed out of the ICM or (2) 
interstellar gas originally in the disk of the galaxy which has been elevated by ram-pressure exerted 
by the fast-moving ICM. The first possibility appears unlikely on the basis that the cooling time 
scale of the ICM (Rangarajan et al. 1995), 

teool ~ 4 X 106(Ticm/10^ i^)^/'(nicM cm-3)-i yr, (1) 

exceeds the Hubble time for the conditions of the Virgo ICM (Tjcm = 2.5 keV = 2 x 10^ K, njcM 
~ 10~^ cm~^; Koyama, Takano, & Tawara 1991; Bohringer et al. 1994). We therefore focus on 
the second scenario. 

We suppose that the extraplanar gas arises from interstellar gas that has been ram-pressure 
stripped from the galactic disk. In this scenario, the ram pressure exerted by the ICM will be 
largest on the receding eastern edge of the disk so this is the most likely source of extraplanar 
gas; this is consistent with the presence of high-|z| gas east of the nucleus. The swept-out gas 
will be heated conductively by the hot ICM; it is therefore important to determine if this gas will 
survive long enough in the ICM environment to be detected optically at the observed large heights 
above the disk. Here, we consider only the HI gas phase because the molecular gas distribution 
in NGC 4388 appears normal (Kenney & Young 1989). Following Cowie & McKee (1977), we 
calculate the evaporation timescale for a cloud embedded in the ICM. In the galaxy halo, classical 
thermal conduction (i.e. no ordered magnetic fields that would reduce the conductivity; see 
Hummel & Saikia 1991) should provide a reasonable estimate for the evaporation time, which is 
given by 

tevap = 1, OOOncRlMcu/W^ K)icM'(lnA/30) yr (2) 
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where ric is the embedded cloud density, Rpc is its radius in parsecs, Tjcm is the temperature of 
the surrounding hot medium, and A is the Coulomb logarithm. For typical Galactic HI clouds (ric 
~ 1 cm~^, R ~ 10 pc) embedded in a hot medium of temperature 2.5 keV (Koyama, Takano, & 
Tawara 1991; Bohringer et al. 1994), the evaporation time is tevap ~ 10'^ yrs. This is much shorter 
than the dynamical timescale tdyn ~ 2 x 10^ years derived in §4.1. 

Shock waves generated during the sweeping may also contribute to heating the interstellar 
clouds. The radial velocity Vr = 1,500 km s~^ of NGC 4388 relative to the ICM would produce 
shock waves in the ICM (where the sound speed cicm ~ 0.1 Tj^j^ km s^^ ~ 500 km s^^) and 
in the interstellar gas of NGC 4388 (chi < 3 km s^^, assuming HI with T < 10^ K). The shock 
velocity in the HI cloud will be Vs ~ (^icm/^c)^" 1^ ~ 15 km s^"*^ (e.g., McKee & Cowie 1975; 
Sgro 1975) and the gas in the cloud will be heated to T ~ 14 w^ ~ 3,200 K in a time t^ ~ '^R/vs 

« 10^ yrs << tdyn- 

Another serious argument against the ram-pressure hypothesis is based on the kinematics of 
the extraplanar gas. The equation of motion for the HI clouds that are stripped by the ICM and 
somehow survive evaporation/heating is 

Mcg = CpicuiVicM - Vcfll (3) 

where g is the cloud acceleration, Mc, T4, and Ic, are the mass, speed, and characteristic size of the 
interstellar clouds, and picM and Viqm are the mass density and speed of the ICM, respectively. 
The value of the constant C depends on the details of the cloud model (e.g., C « 9 for an 
isothermal cloud model; DeYoung & Axford 1969; Longair, Ryle, &: Scheuer 1973; Blandford & 
Konigl 1979a, b). The equation of motion can be simplified if we make the reasonable assumption 
that Vc « ViQM in the present situation and also assume that the size and mass of the clouds are 
constant. (In reality the cloud expands because of the decrease in ram pressure and its mass may 
also decrease because of erosion caused by instabilities). With these assumptions, the distance 
Xacc over which the cloud needs to be accelerated to attain velocity Vc is given by 



V^ ( MrV- 



^^ ) • (4) 



2g y2CpicMV'icM^c 

Using Vc ~ 200 km s^-*^ for the velocity of the extraplanar gas and M^ ~ 20 Mq for the mass of 
typical HI clouds {ric ~ 1 cm^^; Ic ~ 10 pc) entrained in the ICM, we get Xacc of only ^ 100 pc. 

This simple model therefore predicts velocities for the extraplanar material at ~ 2 - 4 kpc 
which considerably exceed those observed. Large projection effects where much of the velocity of 
the entrained gas is in the plane of the sky are unlikely. The location and very large velocity of 
NGC 4388 relative to the Virgo cluster {Vr = 1,500 km s'^ >> dvirgo = 762 km "^ the Virgo 
cluster velocity dispersion; Rangarajan et al. 1995) strongly suggests that NGC 4388 is on a 
radial orbit roughly aligned along our line of sight (otherwise, NGC 4388 would not be bound to 
the Virgo cluster; see also Giraud 1986). The observed tendency for gas at large heights above 
the disk to have velocities close to systemic is opposite to that expected for ISM entrained in the 
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ICM. (To first order, we expect Vc oc |-zp^ in our simple model; eqn. 4). Nor can the brightness - 
velocity correlation in the NE complex and the redshifted velocity of the SW cloud be explained 
in this scenario. We conclude that the extraplanar gas component is not exposed to the putative, 
high velocity ICM. We discuss this surprising result in more detail in §4.3. 



4-2.3. Supernova- Driven Event 

Corbin et al. (1988) have suggested that part of the NE complex (more specifically, the NE 
cloud) may have been elevated above the disk by a recent supernova. This scenario deserves 
serious consideration because evidence for intense star formation is seen throughout the disk of 
NGC 4388 (Fig. 1). Associated with these star-forming complexes are clustered supernovae that 
deposit significant kinetic energy into the surrounding ISM. The kinetic energy from the stellar 
ejecta will thermalize rapidly to create an expanding "bubble" of very hot (T ~ 10^ - 10* K) gas. 
Following Mac Low & McCray (1988) and Mac Low, McCray, & Norman (1989), the condition for 
this bubble to "break out" of the disk is given by 

L> ~ 3 X 10-^2 ^flE/dt) i7fc-2 p-3/2 ^1/2 ^ -^QQ^ (5^ 

where dE/dt is the rate at which energy is injected in the disk, H^pc is the scale height of 
the galactic disk in kpc, uq is the ambient density in cm~'^, and P7 is the ambient pressure, 
P/A;, in units of 10^ K cm~^. In analogy with our analysis of the disk of NGC 3079 (Veilleux, 
Bland-Hawthorn, &: Cecil 1995), we assume that the multi-phase nature of the ISM in NGC 4388 
can be simplified by using the following average disk parameters: Hkp^ = 1.0, no = 1.0, and P-j = 
10^"^. Under these circumstances, break-out occurs if dE/dt ^ 1 x 10'^^ erg s^^. Note that this 
analysis neglects cushioning by the ambient magnetic field (cf. Tomisaka 1990; Ferriere, Mac Low, 
& Zweibel 1991; Slavin & Cox 1992; Mineshige, Shibata, k Shapiro 1993). We have therefore 
probably underestimated the energy injection rate required for break-out. 

Are any of the star-forming regions in the disk of NGC 4388 sufficiently powerful to eject the 
NE cloud out of the disk? Based on the projected location of the NE cloud, the best candidate HII 
region complex for this event lies at i? ~ 1.9 kpc east of the nucleus. This complex has L(IIa) ~ 2 
X 10^^ erg s^^. Using equation (4) of Veilleux et al. (1995) to translate this Ha luminosity into 
a time-averaged rate of energy injection from stellar winds and supernovae, we obtain dE/dt ~ 3 
X 10'^^ erg s~^, barely sufficient for break-out. Note that our analysis has not accounted for the 
kinematic energy in the NE complex itself: £"14111 (NE cloud) = -E'buik(NE cloud) -|- £'turb(NE cloud) 
^ 5 X 10^^ ergs, is quite substantial - the equivalent to the energy released from ;^ 5 Type II SNe. 

Our dynamical analysis suggests that it is difficult to explain the NE cloud as a result of 
star formation. But there are two important assumptions. First, we neglected dust obscuration 
by the edge-on disk when deriving the energy injection rate from the Ha fiux of the eastern HII 
region. Extinction estimates towards the nucleus (Ay ~ 1 - 2 mags; Phillips &: Malin 1982; 
Storchi-Bergmann et al. 1990; Osterbrock et al. 1992; Petitjean & Durret 1993) suggest that 
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the correction may be substantial. Second, it is unclear if the current energy injection rate from 
this complex is relevant to an event that presumably took place 2;(NE)/1/(NE) ~ 10^ years ago. 
Evidently, dynamical arguments alone fail to reject the supernova hypothesis. 

Nonetheless, we consider this hypothesis unlikely for three further reasons: (1) Perhaps 
most damaging is the morphology of the NE complex. It does not resemble the filaments/loops 
produced by intense star formation in the disks of galaxies (e.g., Dettmar 1992; Veilleux et al. 
1995). A bridge of blueshifted gas between the NE cloud and the nucleus is detected in our [O III] 
data. This is strongly suggestive of mass transfer from the nucleus to the NE cloud rather than 
from the disk (see next section). (2) The overall blueshifted velocity field of the NE complex is 
difficult to explain in this scenario. Gas ejected by SN-driven events is expected to have a strong 
velocity component perpendicular to the plane of the disk (e.g., Norman & Ikeuchi 1989). Because 
i(NGC 4388) = -78°, we expect that most of this gas will be traveling away from us, not toward 
us. In addition, the SN-event hypothesis cannot explain the presence of the SW cloud. (3) The 
relatively long orbital period in the disk at i? ~ 2 kpc (2 vr i?/V^rb ~ 10^ years > tdyn) implies that 
morphological (e.g., filaments, bubbles) or kinematical (e.g., non-rotational motion) disturbances 
caused by a supernova should still be visible in the disk. None are detected in our data. 



4-2.4- Nuclear Outflow 

NGC 4388 has an elongated, jet-like radio structure of nonthermal origin which extends ~ 
200 pc south of the nucleus (e.g.. Stone et al. 1988; Carral, Turner, & Ho 1990; Hummel & Saikia 
1991; Kukula et al. 1995; Falcke et al. 1998). This suggests a collimated nuclear outflow. This 
radio morphology resembles that of several other Seyfert galaxies where the ejected radio plasma 
is clearly interacting with the surrounding ISM, producing elongated optical narrow-line regions 
where the gas is compressed and entrained in the outflow, but is probably kept ionized by the 
nuclear continuum (e.g., Whittle et al. 1988). 

Is this phenomenon occurring in NGC 4388? The orientation of the putative nuclear outflow 
in NGC 4388 does not favor strong interaction between radio plasma and ISM in the galaxy disk: 
the nuclear radio source (i.e. the nuclear component and "blob A" in the nomenclature of Hummel 
& Saikia 1991) is elongated along P. A. w 23°, i.e. close to the minor axis of the galaxy disk. A 
diffuse bubble-like structure extends ^ 1.5 kpc above the nucleus due north (P.A. ~ 5°) of the 
galaxy disk. This strongly suggests a nuclear outflow which has burst out of the galaxy plane. The 
existence in Figure 3 of a continuous bridge of approaching (blueshifted) gas between the nucleus 
and the NE complex also supports this picture. The presence of a redshifted cloud south-west of 
the nucleus, i.e. diametrically opposite to the blueshifted NE complex suggests that the outflow is 
bipolar. The mass and kinetic energy involved in the extraplanar gas (cf. §4.1) are similar to those 
involved in galactic-scale outflows in nearby starburst and Seyfert galaxies (e.g., Cecil, Bland, & 
Tully 1990; Heckman, Armus, & Miley 1990; Veilleux et al. 1994; Veflleux k Bland-Hawthorn 
1997; Shopbell & Bland-Hawthorn 1998). Line splitting near the NE and SW clouds may reflect 
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regions of high dissipation of kinetic energy caused by the interaction/deflection of the outflowing 
gas with ambient halo gas (e.g., Galhmore et al. 1996a, 1996b). The blueshifted velocities of the 
extraplanar gas above the north (near) side of the galaxy and redshifted velocities below the south 
(far) side suggests that the radio jet tilts relative to the normal to the disk by at least 90° - |i| Ri 
12°. 

One difficulty with the outflow scenario is the lack of correlation between radio features in 
NGC 4388 and the extraplanar optical gas. Standard plasmon-driven bowshock models predict 
a one-to-one correspondence between the radio and optical structures (e.g. Pedlar, Dyson, & 
Unger 1985; Taylor et al. 1989). While the lack of obvious optical counterpart to radio blob A 
can be attributed to dust obscuration by the disk, the origin of the large offset between the NE 
line-emitting complex in NGC 4388 and the northern radio bubble is harder to explain. The NE 
complex is three times larger than the northern radio bubble, and is misaligned from the axis of 
the radio bubble by ~ 30°. 

However, morphological differences between radio and optical features are common among 
edge-on galaxies with confirmed large-scale galactic outflows (e.g., Colbert et al. 1996a, 1996b, 
1998; Lehnert &: Heckman 1996). In NGC 3079, for instance, most of the radio emission lies 
outside the kpc-scale optical nuclear superbubble (Veilleux et al. 1994), and X-shaped optical 
filaments with no obvious radio counterparts rise more than 4 kpc above the disk plane (Heckman, 
Armus, & Miley 1990; Veilleux et al. 1995). Several processes may explain the lack of correlation 
between radio and optical emission in NGC 4388. The lighter radio plasma is more susceptible to 
external forces than the denser optical filaments. Buoyancy due to density differences between the 
radio cloud and ambient gas in the halo and the galactic gravitational gradient (parallel to the 
minor axis of the disk) have been suggested to explain the bending in the radio structure (Stone et 
al. 1988). Here the net buoyancy force on a cloud of density pc immersed in a region of density po 
is oc pq — Pc- Buoyancy will have no direct effects on the denser optical line-emitting gas, thereby 
kinematically decoupling this gas from the radio plasma. Balancing the buoyancy on the radio 
cloud against the ram pressure from the ambient halo gas and assuming pc « po, Stone et al. 
(1988) derived a terminal velocity for the radio cloud of 10 - 35 km s~^, i.e. much smaller than 
the apparent velocity of the line-emitting gas. Interestingly, the time needed to reach the current 
radio extent at this speed is similar to the dynamical timescale for the line-emitting gas (tdyn ~ 
2 X 10^ yrs). The lighter radio plasma will also be more subject to refractive bending by density 
gradients in the halo gas than the dense line-emitting gas (Henriksen, Vallee, & Bridle 1981; Smith 
& Norman 1981; Fiedler &; Henriksen 1984). A steep vertical density gradient in the ambient 
gas will force the radio plasma to follow a path close to the minor axis of the galaxy, as seen in 
the radio data. Finally, magnetic forces may also act to confine the lighter radio plasma (e.g., 
Wehrle &: Morris 1987) without substantially affecting the kinematics of the more massive optical 
filaments. Current radio data do not constrain the strength and orientation of the magnetic field 
in the extraplanar gas component. 

Assuming the morphological differences between the radio and optical emission are indeed 
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due to a kinematic decoupling between these two gas components, one question still remains: 
is the radio jet in the nucleus of NGC 4388 powerful enough to accelerate the gas in the NE 
complex to its current velocity? To answer this question we assume that the acceleration phase 
of the line-emitting gas occurred near the nucleus when the radio plasma was still coupled to the 
line-emitting gas. The line-emitting clouds, initially at rest, were accelerated downstream by the 
jet. Under these assumptions, the equation of motion of the line-emitting clouds at this early 
epoch is similar to that of clouds entrained in the ICM (eqn. 3). The distance Xacc over which the 
cloud must accelerate to reach velocity Vc is given by equation (4) after replacing picM and Ficm 
by pj and Vj, the mass density and velocity of the jet gas. The value of Xacc can be estimated 
if we use the fact that L^ = eNcPjVH1/2 where L^ is the luminosity of the radio jet, Nc is the 
number of clouds entrained by this jet, and e is the fraction of the total power of the jet which is 
converted into radio emission (Wilson 1981): 
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Here, Mtot = ^c Mc ~ 4 x 10^ Mq is the total mass of line-emitting clouds entrained in the radio 
jet (cf. §4.1). As a measure of the jet radio luminosity, we use the sum of the radio luminosity 
in "blob A" and the "elongated feature" (Hummel & Saikia 1991; L^ ~ 7 x 10^^ erg s^^ for our 
adopted distance of 16.7 Mpc). A reasonable estimate for Xacc is therefore ~ 100 pc, implying that 
the acceleration zone of the line-emitting clouds is well within the radial position of the NE cloud 
(R ~ 2 kpc). Entrainment with the radio jet can therefore produce the observed velocities of 200 
- 300 km s^^. While the previous derivation neglects projection effects and ram pressure by the 
ambient halo gas (both of which will tend to increase the required acceleration distance, Xacc); we 
feel that the kinematic energy in the NE complex can be explained through this nuclear outflow 
scenario. From equation (3), ram-pressure deceleration by the ambient halo gas will be stronger 
for small clouds (large 1^/Mc) than for large clouds, and may therefore explain the larger blueshift 
of the brighter (more massive) clouds in the NE complex. 



4.3. Nature of the Interaction with the Virgo Intracluster Medium 

The disk of NGC 4388 is oriented (i = -78°) such that the extraplanar gas to the north-east 
of the galaxy should feel the full force of the ICM 'wind', while the high-|z| gas to the south 
may be somewhat 'shielded' by the galaxy, lying in the wake of the flow. Simulations of radio 
jet propagation in a moving ICM (e.g., O'Donohue, Eilek, &: Owens 1993; Loken et al. 1995) 
show that the ram pressure from the ICM should have profound effects on the morphology 
and kinematics of the outflowing gas. The lack of evidence for entrainment of the north-east 
extraplanar gas in the Virgo ICM is therefore puzzling. In this section we discuss in more detail 



15 



the nature of the interaction between NGC 4388 and the Virgo ICM, and estimate the effects of 
this interaction on the kinematics of the extraplanar material. 

We first re-examine the rationale for ICM-ISM interaction in NGC 4388. Ram-pressure 
sweeping depends on the local pressure of the intracluster medium, the relative motion of a galaxy 
relative to the intracluster medium and the local surface density within the galaxy. Ram-pressure 
stripping of the gaseous phase of a galaxy will occur if the ram pressure from the ICM exceeds the 
gravitational restoring force in the galaxy (Gunn &: Gott 1972), i.e. 

PicmV± > 2TTGatotcrgas (7) 

where picM is the density of the intracluster medium, V± is the relative velocity of the intracluster 
medium perpendicular to the galaxy plane, o"tot is the galaxy total surface density and cxgas = cthi 
-|- cJHa is the gas surface density. Using V± = 1,200 km s'^, picM = 1-6 x 10~^ cm~^ (derived 
from a general model fit to the X-ray surface brigthness profile; Takano et al. 1989, and references 
therein), cjgas = o"hi = 10^*^ cm~^, and crtot = 157 Mq pc~^ (calculated using a Brandt model 
potential with n = 1.3 at the HI radius where Ugas = 10^*^ cm~^), Cayatte et al. (1994; their 
Table 3) have argued that condition (7) is met beyond R ~ 6 kpc in NGC 4388. But the large 
value of V± used in their derivation and the edge-on orientation of NGC 4388 would imply that 
NGC 4388 has a large velocity transverse to our line of sight, Vt = (y± - Vr cos |i|)/sin \i\ ~ 900 
km s^^. Given the mass of the Virgo cluster, Mvirgo ~ 5 — 9 x IO^^Mq (e.g., Tully & Shaya 1984; 
Giraud 1986), the total velocity of NGC 4388, [Vt^ + V^^l'^ « 1,800 km s~\ would imply that 
NGC 4388 is unbound to the Virgo cluster. We consider this unlikely based on the HI deficiency of 
NGC 4388 and its proximity to the center of the cluster. More likely, NGC 4388 is in a minimum 
potential energy part of its orbit and most of its velocity is along our line of sight. In this case, 
Vj_ ~ T^ cos \i\ ~ 300 km s^^, the ram pressure is l/16th that used by Cayatte et al. (1994), and 
stripping occurs only in the outskirts (R > 10 kpc) of NGC 4388. 

However, as noted by Nulsen (1982), ram pressure stripping is not the only means by which 
the ICM can strip a galaxy. Momentum transfer between the hot ICM and the galaxy ISM via the 
Kelvin-Helmholtz instability can increase stripping rates under some circumstances. While the 
effectiveness of the ram pressure is extremely sensitive to orientation, the rate of mass loss from 
such turbulent viscous stripping is approximately 

Mvisc = ttRhiPicmV ^ -irRuiPicuVr, (8) 

where Rm is the radius of the HI disk and V is the velocity of NGC 4388 relative to the cluster. 
Using the current extent of the HI disk of NGC 4388, Cayatte et al. (1994) estimate that 
-^visc(NGC 4388) ~ 1 - 6 M.q yr^^, implying that it will lose virtually all of its atomic content 
(Mhi = 4 X 10^ Mq for a distance of 16.7 Mpc; Kenney & Young 1989) during a crossing time 
^crossing — 3 X 10^ yrs through the cluster core. This suggests that NGC 4388 is on its first passage 
through the cluster core. 

So far we have assumed that the ISM of NGC 4388 is experiencing the full impact of the 
high-velocity ICM. In reality the ICM 'wind' will be disturbed by the presence of the galaxy. In 
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particular, because the sound speed in the cluster gas is Cg ~ 500 km s^^ and the radial velocity 
of NGC 4388 relative to the Virgo cluster core is Vr ^ 1,500 km s^"*^, the interaction of the ISM 
of NGC 4388 with the cluster gas will be supersonic with Mach number M = V/cs ~ Vr/cg ~ 
3. In this situation, a Mach cone (bow shock) forms near the leading edge of NGC 4388. The 
opening half-angle of the Mach cone is = 2 tan^"*^ (l/M) ~ 40° (note the factor of 2 missing in 
the expression used by Rangarajan et al. 1995). The ICM will be shocked at the Mach cone. The 
Rankine-Hugoniot jump conditions for adiabatic shocks (7 = 5/3) yield a post-shock velocity Vps 
~ ^icm/3 = 500 km s~^, a post-shock density Ups ~ 3 njcM = 3 x 10~^ cm~^, and a post-shock 
temperature Tpg = 3.7 Tjcm = 7.4 x 10^ K. From equation (1), the cooling time for the shocked 
ICM is ~ 3 X 10^^ yrs, thereby justifying our assumption of an adiabatic shock. Given the rather 
wide opening angle of the Mach cone, the extraplanar gas observed in NGC 4388 will be affected 
by the shocked - rather than the pre-shock - ICM (Fig. 5). 

The turbulent viscous stripping rate due to the shocked ICM will be the same as that 
calculated for the pre-shock ICM (because picmVicm = Ppsyps', cf. eqn. 8), but the ram pressure 
on the extraplanar gas will be substantially reduced. Assuming the shocked ICM flow lines roughly 
follow the contours of NGC 4388 and taking K^utflow = 250 km s~^ as the velocity component of 
the outflowing NE cloud along the plane of the galaxy, the ram pressure from the shocked ICM 
on the extraplanar gas will be reduced by a factor pjcu/Pps x [(Vicu - VontHow) / (Vps - Vi^utflow)]^ 
~ 8.3 relative to that from the pre-shock ICM. Under these conditions the acceleration distance, 
Xacc) required to increase the approaching velocity of the gas in the NE cloud from 250 km s~^ to 
350 km s~^ is approximately 



(^final - initial) _ / mVl,, - V,' 
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2<9> \ 2Cpps < (Vps - youtfiow)2 > /2 ; ' 

where V^nai = 350 km s~\ initial = 250 km s~\ and < g > and < {Vps - K)utfiow)^ > are 
time-averaged quantities. Using Mc ~ 20 Mq and Ic ~ 10 pc for typical extraplanar line-emitting 
clouds [i.e., nc ~ 1 cm~^ and Ic ~ Zc(III cloud)], we get Xacc ~ 6 kpc. This is significantly larger 
than the 2-3 kpc separation between the NE cloud and the rest of the NE complex at higher |z|. 
The implication is that the ram pressure exerted by the shocked ICM has only a small effect on 
the kinematics of the outflowing extraplanar material. 

We conclude this section on a more speculative note by suggesting that the impact of the 
shocked ICM on the galaxy and the extraplanar gas component may also be cushioned by a 
gaseous halo in hydrostatic equilibrium around NGC 4388. As we discussed in §4.2.4, a static 
gaseous halo in NGC 4388 would also help to explain possible shock structures in the NE and SW 
clouds and the curvature of the radio jet. Our own Galaxy is known to be surrounded by a hot (^^ 
10^ K) highly ionized corona extending out to several kpc from the Galactic plane (e.g., Spitzer 
1990; Sembach & Savage 1992; Wang 1992). Quasar absorption-line studies (e.g., Lanzetta et al. 
1995) and deep X-ray images of spiral and elliptical galaxies (e.g., Fabbiano 1988; Heckman et al. 
1990; Dahlem et al. 1996) show similar gaseous halos. The main question is whether such halos 
survive long enough in the Virgo cluster environment. Current soft X-ray images of NGC 4388 
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(e.g., Matt et al. 1994; Angelo Antonelli, Matt, & Piro 1997; Iwasawa et al. 1997; Colbert et al. 
1998) show soft X-ray emission of thermal origin, but the exact morphology of the X-ray gas is 
ill-defined because of the strong X-ray background from the Virgo ICM and the distorted halo in 
nearby M86 (Forman et al. 1979; White et al. 1991; Rangarajan et al. 1995). The very existence 
of this halo in M86, a giant E/SO galaxy which is located close to the Virgo cluster core (15.7 Mpc 
according to Jacoby et al. 1990 or 18.4 Mpc according to Tonry et al. 1990) and is moving at a 
velocity ~ -1,300 km s^-*^ relative to the cluster (Vr = -227 km s^"*^; Jacoby et al. 1990), gives us 
some confidence that at least some galaxies in Virgo [albeit those more massive than NGC 4388: 
M(M86) ~ 10^^ M0; e.g., Forman et al. 1979] can keep an X-ray halo while experiencing virtually 
the same ram pressure as that on NGC 4388. 



5. Conclusions 

The complete two-dimensional coverage of our Fabry-Perot spectra on the line-emitting gas 
in NGC 4388 has allowed us to separate kinematically the extraplanar and disk gas components. 
We find that the velocity field of the high-|zj gas is best described by a bipolar outflow rather 
than ram-pressure stripping by the fast-moving ICM in the Virgo cluster. We argue for a Mach 
cone with opening angle ~ 80° to explain these results. Unfortunately, this Mach cone will be 
difficult to detect, even with AXAF and XMM, because it is predicted to be roughly face-on and 
the strong X-ray background from nearby galaxies and the ICM in the Virgo cluster will dominate 
the emission. Nevertheless, high-resolution X-ray spectroscopy would be useful to constrain the 
velocity field (and metallicity) of the X-ray gas and to help establish if the hot gas near NGC 4388 
is in hydrostatic equilibrium in the galaxy potential or instead represents fast-moving, shocked 
ICM. Detailed investigations of bow shocks in galaxy clusters have the potential to help us 
determine the space velocities of galaxy members with respect to the ICM. 

From a theoretical viewpoint there is a need to explore the effects of a fast-moving ICM as it 
interacts with a gas-rich galaxy. Few such models have been published (e.g., Kritsuk 1984; Tosa 
1994), and none considers the multi-phase nature of the ISM and the effects that this interaction 
may have on star formation (e.g., through compression; Kenney & Young 1989; Cayatte et al. 
1994). The ICM 'wind' will also affect the morphology of the nuclear ejecta from starburst 
galaxies and AGN (e.g., wide-angle tailed radio sources in cluster environment; Loken et al. 1995). 
Outflowing gas will stretch in the direction of the ICM motion, thereby increasing the effective 
gaseous cross-section of these galaxies and perhaps contributing significantly to the population of 
C IV absorption lines systems in quasars (e.g., Norman et al. 1996). 
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Table 1: Kinematical Parameters of NGC 4388(") 



Component 


Parameter 


NGC 4388 


Disk 


Vsjs 


2,525 km s^^ 




Kiax (deprojected) 


240 km s-i 




P.A. (major axis) 


90° 




Inclination i 


-78° 


Bar 


P.A. (intrinsic to disk) 


135° 




Inner radius Rh 


1.5 kpc 




Orbit eccentricity R < Rb 


0.3 




Outer radius R^ 


5.1 kpc 



Refer to Veilleux, Bland-Hawthorn, & Cecil (1999) for a detailed discussion of these parameters. 
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Fig. 1.— The distributions of the Une emission in NGC 4388. (top) Ha; (bottom) [O III] A5007. 
The optical fluxes were determined from Gaussian fitting of the emission-line profiles. North is at 
the top and East to the left. The spatial scale, indicated by a horizontal bar at the bottom of the 
image, corresponds to 12", or 1 kpc for the adopted distance of 16.7 Mpc. The optical continuum 
nucleus is indicated in each panel by a cross. The flux scale is logarithmic and is expressed in units 
of 10~^^ erg s~^ cm~^ arcsec"^. 



Fig. 2. — [O III] A5007/Ha flux ratio map of the line-emitting gas in NGC 4388. Same spatial 
scale and orientation as Fig. 1. The optical continuum nucleus is indicated by a cross. Values 
for this ratio range from ~ 0.1 in the disk to ^ 5 in the extraplanar material. Regions with 
[O III] A5007/Hq = 9 indicate [O III] emission clouds where no Ha was detected. 



Fig. 3. — Barycentric velocity fleld of the line-emitting gas in NGC 4388. (top) Ha; (bottom) 
[O HI] A5007. Same spatial scale and orientation as Fig. 1. The optical continuum nucleus is 
indicated in each panel by a cross. The velocites range from about -200 km s^^ to +200 km s^^ 
relative to systemic (= 2,525 km s~^). 



Fig. 4. — Models of the disk velocity field: (top left) a simple axisymmetric kinematic model, 
(top right) a kinematic model with elliptical streaming motion in the center. See text and 
Table 1 for more details on the parameters of the models. The residuals after subtracting the 
axisymmetric model from the observed Ha velocity field are presented in the bottom left panel 
while the bottom right panel shows the residuals after subtracting the non-axisymmetric model. 
The non-axisymmetric model is clearly a better representation of the data. 



Fig. 5. — Geometry of the ICM - NGC 4388 interaction. The values of the parameters indicated 
on this figure and discussed in the text are: i = - 78°, (j) = 40°, 6 > 12°, Ru i = 10 kpc, V[cm 
= 1,500 km s^^ , nicM ~ 10"^ cm^^, Vps = 500 km s^^, rips ~ 3 x 10^^ cm^^. The observer is 
located on the right in the same plane as the figure and at a distance of 16.7 Mpc. 
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